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We study the temperature crossovers seen in the magnetic and transport 

^2 ' properties of cuprates using a nearly antiferromagnetic Fermi liquid model 

00 

O ' (NAFLM). We distinguish between underdoped and overdoped systems on 

<N ' 

the basis of their low frequency magnetic behavior and so classify the opti- 

On 



mally doped cuprates as a special case of the underdoped cuprates. For the 



& ' overdoped cuprates, we find, in agreement with earlier work, mean-field z = 2 

S" 

behavior of the magnetic variables associated with the fact that the damp- 



ing rate of their spin fluctuations is essentially independent of temperature, 
while the resistivity exhibits a crossover from Fermi liquid behavior at low 



^ ■ temperature to linear-in-T above a certain temperature To. We demonstrate 

5-j 



that above To the proximity of the quasiparticle Fermi surface to the magnetic 
Brillouin zone boundary brings about the measured linear in T resistivity. For 
the underdoped cuprates we argue that the sequence of crossovers identified 
by Barzykin and Pines in the low frequency magnetic behavior (from mean 
field z = 2 at high temperatures, T > T cr , to non-universal z = 1 scaling 
behavior at intermediate temperatures, T* < T < T cr , to pseudogap behavior 
below T*) reflects the development in the electronic structure of a precursor 
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to a spin-density-wave state. This development begins at T cr with a thermal 
evolution of the quasiparticle spectral weight which brings about tempera- 
ture dependent spin-damping and ends at T* where the Fermi surface has lost 
pieces near corners of the magnetic Brillouin zone. For < T < T cr the 
resistivity is linear in T because this change in spectral weight does not affect 
the resistivity significantly; below T* vertex corrections act to bring about 
the measured downturn in (p(T) — p(0))/T and approximately quadratic in 
T resistivity for T <C T*. 
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I. INTRODUCTION 



Over the past few years, it has become increasingly clear to the "high-T c " community 
that the mechanism of superconductivity in cuprates is directly related to their unusual 
properties in the normal state, particularly in the underdoped regime. The 63 Cu spin-lattice 
relaxation rate and spin-echo decay rate (Slichter 1994, Barzykin and Pines 1995), uniform 
susceptibility (Johnson 1989), in-plane and c-axis resistivity, p xx and p c (Ong 1990, lye 
1992), all demonstrate temperature dependences which in a wide temperature range are 
different from the predictions of Landau Fermi liquid theory. The remarkable sequence of 
crossovers (from non-universal mean field behavior with a dynamical exponent z = 2 to z = 1 
pseudoscaling behavior to pseudogap behavior) seen in magnetic experiments in the normal 
state of the optimally doped and underdoped cuprates are shown in Figure [[]. Characterizing 
and explaining this behavior, which possesses counterparts in angle-resolved photoemission 
experiments (Shen et al 1996, LaRosa et al 1996, Ding et al 1996), transport measurements 
(Hwang et al 1994), and optical experiments (Puchkov 1996) on the underdoped cuprates, 
is perhaps the major challenge presently facing the high-T c community. 

The effectively non-Landau-liquid behavior at intermediate energy scales has stimulated 
intensive discussions on a possible violation of Fermi-liquid theory in cuprates (Anderson 
1994). There is currently no consensus on whether Fermi-liquid behavior is actually broken 
at T = 0. Some researchers believe that the ground state of underdoped and optimally 
doped cuprates is not a Fermi liquid: in particular, a non-Fermi liquid ground state is a 
point of departure for gauge theories based on spin-charge separation (see, e.g., Lee and 
Nagaosa (1992), Altshuler, Ioffe and Millis (1995)), and for theories in which pairing is 
due to pair hopping between adjacent layers (Chakravarty and Anderson 1994). Another 
conjecture is that there exists an isolated zero-temperature quantum critical point at a 
doping concentration near optimal doping (Varma 1994). On both sides of this point, it is 
proposed that the system behaves as a Fermi-liquid at the lowest T, but possesses a non- 
Fermi-liquid quantum-critical behavior above the crossover temperature, which near optimal 
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doping is assumed to be smaller than T c . 

The authors of the present paper have argued in favor of a third possibility, a nearly 
antiferromagnetic Fermi-liquid model (NAFLM) for cuprates (for a review see, e.g, Pines 
(1995)). In this approach there is no spin-charge separation. One assumes that Fermi-liquid 
behavior is not destroyed by fluctuations at any doping concentration, while the measured 
anomalous (non-Landau-like) spin and charge behavior seen in the normal state properties of 
the cuprates arises from a magnetic interaction between planar quasiparticles which reflects 
the close approach of even an optimally-doped system to antiferromagnetism. Initial support 
for the NAFLM came from the NMR experiments which clearly demonstrate the difference 
between Cu and O spin-lattice relaxation rates and establish the presence of strong tem- 
perature dependent antiferromagnetic correlations (for a review see Slichter (1994)). More 
recently the model has been shown to be also consistent with neutron scattering data which 
show that even at optimal doping, the spin fluctuations are peaked at a wavevector Q which 
is at or near (ir, n) with a magnitude, Xq ^* Xo, an d a half width, l/£, which at low temper- 
atures is considerably smaller than the inverse lattice spacing (Mason et al 1996, Bourges 
et al 1996, Zha, Barzykin and Pines 1996). Moreover, recent experiments have shown that 
even at optimal doping there still exist propagating spin-waves at energies comparable to 
the exchange integral (Hayden et al 1995). 

In this paper, we consider the extent to which the NAFLM can explain the temperature 
crossovers in magnetic and resistivity measurements in both overdoped and underdoped 
cuprates. We will show that the NAFLM yields a sequence of crossovers and changes in 
the uniform susceptibility, NMR relaxation rates and resistivity which are consistent with 
experiments. In particular, we will show how for the underdoped cuprates the Fermi surface 
evolution which accompanies the development of a precursor to a spin-density-wave state 
gives rise to z — 1 scaling at intermediate temepratures, and the pseudogap behavior at the 
lowest temperatures. 

The paper is organized as follows: in the next section we review the NAFLM description 
of cuprates and its relation to the underlying microscopic models with fermion-fermion 
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interaction. In Sees. |T|, [TV] and [V] we discuss the physics of overdoped, underdoped and 
optimally doped cuprates, respectively. Finally, in Sec. |VJ we summarize our conclusions. 



II. THE NEARLY ANTIFERROMAGNETIC FERMI LIQUID MODEL 

The canonical NAFLM model is in some respects a two-fluid model for cuprates: it is 
equivalent to assuming that there independently exist fermionic carriers and localized spins 
whose susceptibility is determined from fits to NMR experiments and is an input parameter 
in the theory. Fermions and spins are coupled by 

Him = #q c k+q,a^a,/3Ck,/3 • S q (1) 

q,k,a,/3 

where g q is the momentum-dependent coupling constant, and (7j are the Pauli matrices. The 
momentum dependence of the coupling constant is not relevant for our considerations, and 
for simplicity we will neglect it throughout the paper. 

To second order in the coupling, the spin-fermion interaction gives rise to an effective 
pairing interaction between planar quasiparticles 

V eff (q,uj) = g 2 x((l,uj) (2) 

where x(q, w ) is the susceptibility of localized spins. NMR and neutron scattering exper- 
iments clearly indicate that near optimal doping the spin fluctuations are overdamped at 
low energies with a dynamical structure factor peaked at a wave vector Q which is close to 
(71", 7r) and symmetry related points. From general considerations, one can then write near 
q = Q and uj = (Millis, Monien and Pines 1990, Barzykin et al 1993, Monthoux and Pines, 
1994b): 

x(<l> u ) = TT~( XQ (3) 

1 + (q - Q) 2 £ 2 - iuj/u s f - uj% 2 /cf 



where xq = a£ 2 , £ is the correlation length, c sw is the spin-wave velocity, and u s f = c 2 sw /2^ 2 ^ 
where 7 is a damping rate. The fits to the NMR data show that the parameters, £, xq an d 
u s f depend on temperature, while the scale factor a and the spin-wave velocity c sw are 



virtually independent of T in the temperature range of NMR experiments (Barzykin and 
Pines 1995). 

It is also essential that the fermions are not assumed to be free particles at g — 0, and their 
Green's function is well defined only near the Fermi surface where G(k,u) = Z^/(uj — e^). 
The dispersion, e^, is another input parameter in the theory. It is generally assumed to have 
the same form as in the tight-binding model with hopping between nearest and next-nearest 
neighbors on a square lattice: 

6k = — 2t(cos k x + cos k y ) — 4t' cos k x cos k y (4) 

This form of e& (with t' < 0) is chosen to be consistent with the shape of the Fermi surface 
measured in photoemission experiments at around optimal doping (Si et al 1993). Explicit 
calculations show (Monthoux and Pines 1993) that at these doping levels, the interaction 
with spin fluctuations does not modify substantially the shape of the Fermi surface as long 
as the coupling constant remains roughly smaller than the quasiparticle damping rate. 

In principle, the correlation length in (^) should be large enough compared to the inter- 
atomic spacing to ensure that the susceptibility is peaked at Q. In practice (Pines 1995), 
it turns out that a correlation length of the order of a lattice spacing already yields an ap- 
preciable peak in x(q, to). Furthermore, the analytical expansion in (q — Q) 2 and u is valid 
only outside the fluctuation region for a magnetic phase transition. This region is however 
rather narrow and most experiments are performed outside it. A more essential point is 
that x(q> u ) m Eq. (|3|) does not satisfy the sum rule for the local susceptibility, as the 3D 
integral / d 2 qdux{<i, uj) diverges at the upper limit. One way to obtain convergence is to 
include higher-order terms in the expansion over momentum. Another way, which we will 
adopt here, is to introduce a sharp cutoff, C, in the integration over c sw (q — Q). Due to 
the thermal dependence of £ and u s f in Eq. (|J), the cutoff parameter C also depends on 
temperature. On general grounds, the larger the correlation length is, the smaller is the 
cutoff scale. However, C does not vanish when the correlation length diverges. 

Although the phenomenological NAFLM considers localized spins as an independent 
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degree of freedom, Eqs. (||) and (Q) can be derived from microscopic considerations departing 
from, e.g., the one-band Hubbard model which contains only fermionic degrees of freedom 
(Shraiman and Siggia 1988, Schrieffer, Wen and Zhang 1989, Chubukov and Frenkel 1992, 
Bulut, Scalapino and White 1993, Scalapino 1994, Kampf 1995, Sachdev, Chubukov and 
Sokol 1995). In this approach, spin fluctuations appear as a collective mode of fermions. To 
obtain the coupling between two fermions and one spin fluctuation, as in Eq. (|l|), one has to 
dress the original four-fermion Hubbard interaction term by summing the RPA series in the 
particle-hole channel. As a result of the summation, the product of the two fermionic Green's 
functions is replaced by a spin susceptibility whose poles correspond to spin fluctuation 
modes. The susceptibility thus obtained takes the general form 

in which x(q, u) is the irreducible particle-hole susceptibility, and g (the same as in Eq. 
(|I|)) is a coupling constant for the pair interaction between fermions. For example, for 
on-site Hubbard interaction, g = U . The momentum and frequency dependence of the 
spin susceptibility obtained in the Hubbard-based approach is generally consistent with 
the phenomenological predictions except that at any finite doping level the peak in the 
susceptibility is located at an incommensurate momentum Q determined by the maximum 
in x(q, 0). (see, e.g. Griiner 1994). It is essential however that for fermionic dispersion as in 
Eq. (|4]), the enhancement of x(q, ^) comes solely from near resonance in the denominator in 
Eq. (|5|), while the irreducible particle-hole susceptibility does not contain any information 
about strong magnetic fluctuations, and near the peak at q = Q has a conventional Fermi- 
liquid form 

m,<*) = z a ~ e lf , (6) 

1 - IUJ/Yq 

where £ is of the order of the lattice spacing, and Tq is an energy comparable to a fermionic 
bandwidth. We also assumed that there exist points on the Fermi surface which can be 
connected by Q; otherwise, Tq would be infinite. It then follows from Eq. (||D that 



sw 
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We see therefore that the key feature which makes the physics in the nearly-antiferromagnetic 
Fermi liquid different from that in a conventional Fermi liquid is the presence of an energy 
scale u s f oc £~ 2 , which vanishes when magnetic correlation length becomes infinite, and is 
obviously much smaller than the fermionic bandwidth as long as the susceptibility is strongly 
enhanced near Q. 

There are, however, several subtleties with the microscopic derivation of Eq. (|3[). First, 
in a RPA calculation, the fermions which contribute to the spin susceptibility are generally 
assumed to behave as free particles everywhere in momentum space, i.e., one completely 
neglects the incoherent part of the quasiparticle Green's function. For the calculations of 
the imaginary part of x{% u), this assumption is justified as the integration over fermionic 
momentum is confined to a region where both fermions in the bubble have momenta near 
the Fermi surface. For these fermions, the incoherent part of the fermionic Green's function 
is irrelevant. The special points in momentum space for which both k and k + Q are on the 
Fermi surface are usually referred to as hot spots (Hlubina and Rice 1995). The real part 
of Xi on the other hand, comes from an integration over regions in momentum space which 
are far from the Fermi surface (see, e.g., Rickayzen 1980). In these regions, the incoherent 
part of G(k,u) cannot be neglected as demonstrated in photoemission experiments (Wells 
et al 1995, Campuzano et al 1994). It is therefore likely that the actual position of the peak 
in Rex{q, is different from the RPA result and may well depend not only on doping but 
also on temperature. 

Second, even though the computation of Imx{q, uj) involves only coherent parts of the 
quasiparticle Green functions, the RPA approach neglects possible strong vertex corrections 
to the polarization bubble. Explicit calculations show that these corrections are the strongest 
for the hot spots (Chubukov 1995, Altshuler, Ioffe and Millis 1995, Amin and Stamp 1996, 
Monthoux 1996). At small coupling they are obviously small to the extent that g / y/juisf < 1. 
At somewhat larger couplings, g > ^fuJJJ, the relative vertex corrections which come from 
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the integration near the Fermi surface scale as [g 2 u s fXQ/v F ] log(C/o; s /), where vf is the 
Fermi velocity (we set the interatomic spacing a = 1). Now, if the damping term is com- 
puted self-consistently, then uj s f by itself scales as g~ 2 , and the relative vertex correction 
depends on g only logarithmically. More explicitly, we have g e ^ = g[l+2/3 log(C/u) s f)} where 
(3 depends only on the shape of the Fermi surface near the hot spots, and is about 1/16 for 
the experimentally measured Fermi surface at optimal doping (Altshuler, Ioffe and Millis 
1995, Chubukov 1995). Notice that vertex corrections increase the spin-fermion interaction 
(g e ff > gr) and hence act in favor of magnetically-induced, Eliashberg-type d— wave super- 
conductivity. We see that f3 is small numerically, so that if the logarithm is not large, the 
relative vertex correction can be neglected. In this case, the RPA analysis is approximately 
valid, and the imaginary part of the full susceptibility comes solely from the imaginary part 
of the particle-hole bubble. This is what we believe happens in the overdoped regime, which 
we associate with moderate values of g ~ y/juj s f in Eq. (p. However, we will argue below 
that the underdoped cuprates are described by a spin-fermion model with a somewhat larger 
ratio gj ^J^uj s f, where new physics associated with Fermi surface evolution appears. If we 
formally keep using the self-consistent approach, we find that at the onset of the Fermi sur- 
face evolution, the logarithmic term overshadows the smallness of f3, and the relative vertex 
correction is not small. In particular, for g = g£) defined below, 2/3 log C /u> s f ~ 7r/4. In 
this situation, the validity of the self-consistent RPA calculations of the damping term also 
becomes questionable. 

From the above considerations we see that, although some properties of the NAFLM 
can be derived in Hubbard-based calculations, one has to be careful about using an RPA 
formalism in calculating the dynamical spin susceptibility. We therefore advocate a semi- 
phenomenological approach in which it is assumed that the spin susceptibility has the form of 
Eq. (0), with both £ and 7 taken from fits to the experimental data. We emphasize however 
that from a physical perspective, only £ should be considered as an independent input 
parameter for Eq. ([I]). The damping of spin excitations due to the interaction with fermions 
is fully described by a spin-fermion model, and we will use the phenomenological form for 
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u s f only because we are currently unable to self-consistently compute the spin damping at 
g ^> c s „,£ _1 . We will, however, be able to estimate the value of the spin damping at very 
large couplings where the electronic structure develops the precursor of a spin- density- wave 
form. 

We now consider the extent to which the physical properties of the normal state can be 
explained by the NAFLM. 

III. OVERDOPED CUPRATES 

We first discuss the overdoped cuprates. We define overdoped systems as those whose 
magnetic behavior is such that in the normal state their uniform magnetic susceptibility, Xo ; 
is either independent of T or even increases with decreasing temperature, while the product 
63 TiT and the square of the 63 Cu spin-echo decay rate ( 63 T 2 g) 2 are linear in T. Examples 
of overdamped systems are L^-zSr^CuC^ for x > 0.2 and Bi 2212 and Tl 2201 for the 
appropriate oxygen content. The measured resistivity of overdoped cuprates is linear in T 
above some particular temperature T = T and has a Fermi liquid form p ~ T 2 at T <C T . 

Our theoretical interpretation of the experiments on overdoped cuprates is based on the 
conjecture, first put forward by Barzykin and Pines (1995), that in these systems 7 and c sw 
are independent of temperature, while £ < 2 for all T > T c . This in turn implies that for 
overdoped cuprates 

uj s f£, 2 — const. (8) 

Since the correlation length £ (which, we remind the reader, is an input parameter for 
NAFLM) does not exceed a few interatomic spacings, its temperature dependence is governed 
by fluctuations at lattice scales and is likely to be material dependent. NMR experiments 
tell us that for all compounds studied £~ 2 scales linearly with T, i.e., £ -2 = A + BT, where 
A and B are (material dependent) constants. More specifically, fits to the NMR data using 
Eq. (£J) and the standard Shastry-Mila-Rice Hamiltonian for hyperfine interactions (Shastry 
1989, Mila and Rice 1989), then yield (Thelen and Pines 1994): 
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in agreement with the experimental results. The form of the susceptibility with 7 and c sw 
independent of temperature and £~ 2 = A + BT is formally equivalent to that found in the 
z = 2 quantum-critical regime (Sachdev, Chubukov and Sokol 1995). Clearly then, the T 
dependence of all observables associated with the spin susceptibility near (71", it) should be 
the same as for z = 2 scaling (Millis, Monien and Pines, 1990). In particular, the measured 
ratio 63 TiT/ 63 T| G oc w s ^ 2 is independent of temperature. 

However, the analogy with z = 2 scaling is only a formal one. Although it works for 
observables, such as 63 T! and T 2 g, which probe the susceptibility near (71", 7r), the behavior 
of the uniform susceptibility does not follow the scaling prediction, which is based on the 
assumption that spin fluctuations near q = and q = Q are related. These are not related 
for the overdoped cuprates because their physics is dominated by the short wavelength 
fluctuations. Indeed, above T cr , Xo(T) decreases with increasing T (NMR studies show that 
T cr roughly corresponds to the temperature where Xo(T) is a maximum). On the other 
hand, 1/N quantum-critical calculations yield Xo(T) linearly increasing with T in the z = 2 
regime (Sachdev, Chubukov and Sokol 1995) [] As microscopic self-consistent (Eliashberg) 
calculations show (Monthoux and Pines 1993, 1994a), the decrease of susceptibility with 
increasing temperature is a natural consequence of lifetime effects, which act to reduce 
Xo(T), and which play an increasingly important role as the temperature increases. 



1 There is some uncertainty in the scaling prediction for Xo(T) because d = 2 is an upper critical 
dimension for the T = transition. Ioffe and Millis argued that the linear term in susceptibility 
could be of either sign (Ioffe and Millis 1995). The sigma- model based calculations yield a positive 
linear in T term (Sachdev, Chubukov and Sokol 1995). 
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The behavior of the resistivity p is also nonuniversal, and, in principle, depends on 
the details of the electronic band structure. However, as we now show, in an argument 
which we shall see applies to the optimally doped cuprates as well as overdoped cuprates, 
p(T) is in fact proportional to T in a wide temperature range due to the proximity of the 
quasiparticle Fermi surface to the magnetic Brillouin zone boundary. The reason is the 
following: according to the Drude formula, which is applicable to the overdoped cuprates 
because vertex corrections are small, the conductivity is proportional to the quasiparticle 
relaxation time averaged over the Fermi surface. The inverse relaxation time - the relaxation 
rate - is given by the imaginary part of the quasiparticle self-energy. As discussed above, 
the overdoped regime is associated with small to moderate values of the coupling constant. 
In this situation, second-order perturbation theory is likely to be valid. To second order in 
g we find by performing a summation over frequencies: 



where n(e) and /(e) are Bose and Fermi distribution functions, respectively. The integration 
over k' can be split into an integration over and an integration over surfaces of equal en- 
ergy. For k located at the Fermi surface, the first integral can be evaluated analytically. The 
exact result is rather cumbersome, but to a good numerical accuracy it can be approximated 

as 



where Ukk' = <^s/(l + (k — k' — Q) 2 £ 2 ), and k and k! are two points at the Fermi surface. 
The same expression is indeed obtained in an ordinary Fermi liquid. In the latter case u s f, 
and hence u^y are of the order of the Fermi energy which is large compared to T. One 
then immediately obtains l/r^ oc T 2 for all k, i.e., the resistivity is proportional to T 2 as 
expected. In cuprates, however, u s f is much smaller than the Fermi energy; the dominant 
contribution to the integral over k! in ( |T2"D then comes from the regions where the distance 
between k' (at the Fermi surface) and k — Q (at the "shadow" Fermi surface) is minimal. 




(11) 




(12) 
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Clearly, this minimal distance is zero at hot spots. Expanding near each of the hot spots, 
and performing the integration over k', we obtain 

T k = 2C 7^ l 1 ^) 

where 



u; a/ (l + (Ah)*?) + + u, f (l + {Akye)) , (14) 

Here Ak is the displacement of k from a nearby hot spot, and v F s is the Fermi velocity at 
the hot spots. It is clear from Eq. (|14]) that r k is the largest in the regions of the Fermi 
surface where Ak is greater than the inverse correlation length, i.e., (Ak) 2 !; 2 ^> 1. These 
regions provide the dominant contribution to the conductivity, while the regions near hot 
spots contribute very little. 

On averaging Qk over the Fermi surface, one finds that < Qk > is independent of T for 
T < T , and < Qk >~ T^/uJJf^ for T > T , where the crossover temperature, T , is given by 

T w u sf (Ak max ) 2 e/2n. (15) 

The extra factor of \ in T arises for numerical reasons. In general, one might expect that 
T is of the order of the hopping integral, i.e., it is much larger than T over the whole 
experimentally probed range of temperatures. If this was the case, then the above reasoning 
would imply that the resistivity, p ~ T 2 £^/uj s f/ < Q k >, would be roughly quadratic in 
T. However, photoemission experiments performed at or near optimal doping show that a 
substantial portion of the Fermi surface is located near the magnetic Brillouin zone boundary. 
It is not unreasonable to assume that this will be the case for the overdoped materials as 
well. In this case To turns out to be much smaller than the hopping integral. For example, 
in YBa2Cu307 on finds T ~ 15meV-C t ~ 250meV. In this situation < Qk >~ Ty/uJJp^ 
starting already from relatively low temperatures and hence p(T) ~ T. Our numerical 
analysis performed for the parameters of the NAFLM chosen to fit NMR experiments in 
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YBa 2 Cu307 (Monthoux and Pines 1994a; Stojkovic 1996) shows that the resistivity is indeed 
roughly proportional to T over the temperature range T > T ~ 15meV. Hlubina and Rice 
(1995) obtained To ~ 20meV with the same values of spin-fluctuation parameters, but 
somewhat different band parameters; hence they found p ~ T 2 up to higher temperatures. 
We emphasize once again that the linear behavior of the in-plane resistivity is due to the fact 
that uj s f is the smallest parameter in the problem, even though it scales as T: the typical 
u s f(Ak) 2 !; 2 are smaller than 2nT, but much larger than u s f. If ou s f oc T was comparable 
to ttT above T cr , then the resistivity would be proportional to T 1 / 2 . Alternatively, if u s f 
is large and nearly temperature independent, as is the case in the more heavily overdoped 
cuprates, then p would be proportional to T 2 , as seen experimentally. Notice also that, since 
< Qk >oc T^JuJJj^ above T , the resistivity does not depend on the spin-damping and hence 
on co> s /£ 2 . We make use of this fact below when we discuss the behavior of resistivity at 
smaller doping concentrations where 7 acquires a strong T dependence. 

It is also instructive to discuss the temperature dependence of the average scattering 
rate, <C l/r^ ^> for T > T . Under the same conditions as above, Qk oc T(Ak), and the 
integration over Ak yields < l/r k >oc Tlog((AA; ma:E ) 2 /u; s /). However, the T dependence 
coming from the logarithm is rather weak and to a good accuracy < 1/r^ >oc T. Therefore, 
we see that at T > T cr , inverting the average of r k and the averaging of l/r^ yields a similar 
linear in T dependence of the resistivity. 

IV. UNDERDOPED CUPRATES 

We now turn to the key issue of the paper, the behavior of the underdoped materi- 
als. These may easily be distinguished from overdoped materials on the basis of their low 
frequency magnetic behavior, depicted in Fig. [L] (Barzykin and Pines 1995). Examples of 
underdoped cuprates include YB^CuziOs, YB^CusOe+x, for x < 0.93, L^-xSr^CuC^ for 
x < 0.2, Bi 2 Sr 2 CaCu208-<5 (Bi 2212) compounds with adjusted oxygen partial pressure dur- 
ing annealing and Hg 2223. As seen in Fig. [I], underdoped cuprates exhibit two different 
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crossovers in the normal state at the temperatures T cr and T* < T cr . Above T cr their behav- 
ior is similar to what has been observed for overdoped systems: it is nonuniversal, but the 
resistivity is approximately linear in T, and the ratio 63 TiT/ 63 T 2 2 G is independent of T, as 
shown for YI^CXuOs in Fig. ||a. The linear behavior of the resistivity persists below T cr 
with only a subtle change in slope; however, the magnetic properties change rather drasti- 
cally below T cr . As shown in Fig. [|b for YE^C^Os, 63 Ti becomes almost independent of T, 
while T2G becomes proportional to T in such a way that the ratio 63 TiT/ 63 T2G is independent 
of T. Using Eqs. @ and fllCf ), we find that in this T range the inverse correlation length 
and uj s f both increase linearly with increasing T, in such a way that ou s f£ — const, implying 
z = 1 behavior. The uniform susceptibility is also proportional to T below T cr . Finally, 
at even lower temperatures, T < T*, the system enters into a regime where the correlation 
length becomes independent of T, m T{T oc u s f increases as T decreases further towards T c 
( 63 TiT displays a minimum at T ~ T#), and both the uniform susceptibility and l/ 63 TiT 
fall off sharply with decreasing T. This regime has been called a "pseudogap regime" as 
the behavior of, e.g., Xo(T), is, at first sight, reminiscent of the behavior of systems which 
display a true quasiparticle energy gap, such as a superconductor or an ordered antiferro- 
magnet. However, the curvature of the fall-off in Xo(T) is opposite to that found below T c 
in conventional superconductors (Slichter 1994). 

We see therefore that the new features of the underdoped cuprates, as illustrated in 
Figure [I], are (z) the crossover from z = 2 to z = 1 behavior at T cr , and (ii) the crossover to 
pseudogap behavior at an even lower T*. 

Before presenting our scenario for these crossovers, we briefly review the results of the 
o— model based studies of the low-temperature crossovers between different universal scaling 
regimes in disordered antiferromagnets (Chakravarty, Halperin and Nelson 1989, Chubukov, 
Sachdev and Ye 1994, Sachdev, Chubukov and Sokol 1995). In the disordered state at T = 0, 
the antiferromagnets possess a gap, A, in the excitation spectrum. If the damping term at 
T = is much smaller than A, then, as T increases, the system experiences a crossover from 
a quantum-disordered regime with an exponential (e~ A / T ) behavior of observables, to the 
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z — 1 quantum-critical regime. If, on the contrary, at T = 0, 7 ^> A, then the sequence of 
crossovers with increasing T is from a quantum-disordered regime with a power-law behavior 
of observables to the z = 2 scaling regime with overdamped spin fluctuations, and finally, 
at even larger T, to the z — 1 scaling regime in which typical magnon frequencies and the 
damping term both scale linearly with T. 

We see that for all ratios of 7/ A the prediction of the sigma-model based studies is 
that z — 1 behavior always occurs at higher temperatures than z = 2 behavior. This is 
simply due to the fact that z — 1 scaling requires that typical frequencies should be larger 
than the T = value of the damping term. Experimentally, however, the situation is the 
opposite: in underdoped cuprates, one observes purely relaxational behavior at high enough 
temperatures, and z = 1 scaling at lower temperatures (Barzykin 1996). The crossover at 
low T from the quantum-disordered behavior to the z = 1 quantum-critical behavior is more 
consistent with the observations around T*, but again, the opposite curvature of the fall-off 
in Xo(T) compared to the sigma-model prediction indicates that the crossover at T* involves 
not only localized spins, but also electronic degrees of freedom. 

Monthoux and Pines (1994b) proposed that the reversed behavior observed in under- 
doped cuprates below T cr is due to the strong temperature variation of the damping term 7 
in the dynamical spin susceptibility, which in turn gives rise to the anomalous temperature 
dependence of u> s f. This temperature variation is neglected in the a-model approaches which 
assume that throughout the whole temperature range of interest, the input parameters in 
the dynamical susceptibility retain the same values as at T = 0. Below we relate the temper- 
ature variation of the damping rate to the formation of a precursor to a spin- density- wave 
state which in turn causes the evolution of the quasiparticle Fermi surface with temperature 
and doping concentration. We argue that the two crossover scales T cr and T*, observed in 
underdoped cuprates indicate the onset and the end point of this Fermi surface evolution, 
respectively. At T cr , the onset of the crossover, the quasiparticle residue along the Fermi 
surface develops a minimum at hot spots; at T*, the end point of the crossover, the system 
actually begins to lose pieces of the Fermi surface. In other words, we argue that there exists 
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only one extended crossover in system behavior in which the electronic structure gradually 
develops the features of a precursor to a spin-density-wave state. In this scenario, z — 1 
scaling in the temperature range T* < T < T cr is just an intermediate asymptotic in the 
extended crossover region, rather than an extension of the z — 1 quantum-critical scaling 
observed in the intermediate T range right at half-filling. 

Fermi-surface evolution with doping at fixed T has been observed in photoemission ex- 
periments on YBCO (Liu et al 1992) and more recently in experiments on Bi-2212 (Marshall 
et al 1996, LaRosa et al 1996). These experiments demonstrated that while near optimal 
doping the hole Fermi surface is large and encloses an area consistent (to within the accuracy 
of the measurement) with the Luttinger theorem, the measured Fermi surface in underdoped 
materials looses pieces near (n, 0) and symmetry-related points. Recent experiments by Ding 
et al (1996) have shown that the same effect occurs at a fixed doping as the temperature is 
lowered. The transformation of spectral weight from the low frequency part of the spectrum 
to higher energies upon approaching half-filling has been observed by Timusk and his col- 
laborators in optical experiments on the planar conductivity a(cu,T) (Puchkov et al 1996); 
the same experiments also demonstrated that there is no spectral weight transformation at 
optimal (and larger) doping. 

Recently, one of us has considered Fermi-surface evolution with increasing spin-fermion 
coupling constant in the NAFLM at T = (Chubukov, Morr and Shakhnovich 1996). 
It was found that as g increases, the quasiparticle residue, Z near hot spots decreases. 
The decrease of Z becomes appreciable when the coupling constant g exceeds a critical 
value g£) ~ vp/(uj s fXQ l°g(C/ w s/)) 1//2 - Notice that g£) vanishes logarithmically when the 
correlation length becomes infinite. As g increases even further, the quasiparticle Fermi 
surface undergoes a substantial evolution in the process of which parts of the Fermi surface 
near the corners of the magnetic Brillouin zone (where the hot spots are located when g ^ 0) 
move away and, simultaneously, there appear two distinct peaks in the density of states - the 
precursors of the valence and conduction bands. This evolution of the Fermi surface occurs 
in the range of g values of the order of the upper cutoff in the staggered spin susceptibility, 
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g ~ g$ ~ C > g£). It was further argued that as the system approaches half-filling, the 
effective coupling constant g increases, while both gj£) and gf) decrease, so that in varying 
the ratio g/g cr one in fact varies the doping concentration. 

Suppose now we fix the doping and vary the temperature. As the temperature decreases, 
the inverse correlation length and the damping rate decrease, and, hence, both critical values 
of g go down (g£) goes down chiefly because of the decrease of damping, while gj® goes down 
chiefly because of the decrease of As a result, even if the effective coupling constant 

weakly depends on T, the ratios, g/g 1 ^) and gj gf) still increase with decreasing T. This in 
turn implies that, as the temperature is lowered, one should observe crossovers in the system 
behavior analogous to those described above. These crossovers will occur at temperatures 
at which g£' 2 \T) become equal to a given g. It seems natural to associate the temperature 
at which g^)(T cr ) = g with T cr , and the temperature at which gV)(T*) = g with T*. 

We now consider the conditions under which the system possesses z — 1 scaling behavior 
below T cr . Scaling with z — 1 requires that typical frequencies and momenta be of the order 
of T. In the NAFLM, this will be the case if both and ou s f scale linearly with T, which 
in turn implies that Tq^ and 7 should exhibit the same linear dependence on T as does 
(Chubukov, Sachdev and Ye 1994, Monthoux and Pines 1994b). Thus one should have 
between T cr and T*, 

7~fQ i ~r 1 = a+ Kr ( 16 ) 

while, as we have seen, above T cr , 7 and Tq 1 are independent of T. As we discussed above, the 
temperature-dependent inverse correlation length should be considered as an input parame- 
ter in the NAFLM because if one attempts to compute the real part of the spin susceptibility 
(from which £ is inferred) in RPA-type calculations, one would need to know the exact form 
of the fermionic Green's function far away from the Fermi surface. Hence, the linear in T 
dependence of £ _1 below T cr follows from the fit to the NMR data on 63 Ti and T 2 (Barzykin 
and Pines 1995). On the contrary, the linear T dependence of 7 has to be obtained within 
the NAFLM approach, although at the moment we do not have a clear recipe for calculating 
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it. 

It is essential however not to confuse the intermediate z — 1 scaling regime in underdoped 
cuprates with the z — 1 quantum-critical scaling in pure antiferromagnets. In the latter case, 
the damping of spin excitations is due to the interaction between spin fluctuations, and one 
can argue quite generally that it should be linear in T when oc T (Chubukov, Sachdev 
and Ye 1994). However, the estimates of Barzykin and Pines (1995) show that at doping 
concentrations where experiments on underdoped cuprates have been performed, the purely 
magnetic damping is too small to account for the experimental data. In other words, in the 
z = 1 regime in underdoped cuprates, the damping primarily comes from the interactions 
with fermions (this is why we describe this regime as "pseudoscaling" ) . At the same time 
we expect that below T cr , lattice corrections are not large. Then, once the z — 1 form of the 
dynamical susceptibility for q near Q is obtained, one can use the universal scaling forms 
obtained for z — 1 quantum-critical scaling for both uniform and staggered susceptibility. 
In particular, this implies that the static uniform susceptibility should be linear in T, in 
agreement with experiment. This, we recall, is different from the z = 2 regime at higher 
temperatures, where the physics is dominated by short wavelength fluctuations and one 
could only formally use the predictions of the z = 2 quantum-critical scaling to obtain the 
T dependence of observables associated with the dynamical susceptibility near Q; however, 
the temperature dependence of the uniform susceptibility was very different from the scaling 
prediction. 

We note in passing that the fact that spin damping is chiefly due to the interaction with 
fermions may resolve an apparent discrepancy between the values of the spin-wave velocity 
in YBCI2CU3OQS3 extracted from a linear fit to the uniform susceptibility and from the 
ratio of 63 T\T ' /T 2G (Millis and Monien 1994, Chubukov, Sachdev and Sokol 1994). Indeed, 
the slope of the uniform susceptibility xo oc T is finite in the absence of spin damping. 
It does not change much with doping compared to its value at half-filling, and therefore 
is likely to only weakly depend on 7. Accordingly, the fit for Xo(T) yields a c sw which is 
slightly larger than that found at half-filling, which agrees with analytical calculations at 
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small doping (Chubukov and Frenkel 1992, Schulz and Zhou 1995). At the same time, the 
NMR relaxation rate is proportional to the damping, and for this quantity it is essential to 
know where the damping comes from. A fit to NMR data using the a-model results yields 
a spin-wave velocity which is a few times smaller than at half-filling. If instead, we extract 
c sw from the NMR data using Fermi-liquid damping, we obtain a spin-wave velocity larger 
by a factor 7fl/7s/ which better agrees with the velocity extracted from the susceptibility 
data. 

We consider next the behavior of the resistivity. Experimentally, it is linear in T both 
above and below T cr with about the same slope. On the theoretical side, we have seen above 
that the resistivity is linear in T as long as over substantial portions of the Fermi surface, 
T is larger than T . This result holds approximately even when 7 acquires a substantial 
T dependence, provided the vertex corrections are not too large. More importantly, the 
decrease of the spectral weight near the hot spots has little effect on the conductivity, 
simply because the dominant contributions to conductivity come from the pieces of the 
Fermi surface which are relatively far from hot spots. A careful inspection of the measured 
shape of the Fermi surface in YBCO compounds (Gofron et al 1994, Campuzano et al 1994) 
and the fits to NMR measurements (Barzykin and Pines 1995), shows that in the underdoped 
cuprates To is also substantially smaller than T cr . In the intermediate, z — 1 scaling regime, 
7 acquires a temperature dependence, and T becomes a function of T, scaling as T ~ T -1 . 
However, the slope is rather small and everywhere in the z = 1 regime T remains smaller 
than T, i.e., the resistivity remains linear in T, in good agreement with the experimental 
results in, e.g., YBa 2 Cu 4 8 (Bucher et al 1993, see Fig. |). The two temperature scales 
(T and To(T)) become comparable to each other at a temperature which in YBa2Cu40s 
is rather close to T* (T* is roughly 220K and T (T^) ~ 250K). We emphasize, however, 
that the behavior of resistivity in our model is nontrivial and depends on the shape of the 
quasiparticle band structure. Hence our argument, based on YBa 2 Cu40s, that the crossover 
in resistivity occurs near T* need not necessarily apply to all high-T c compounds. 

We now turn to the crossover at T = T*. Barzykin and Pines (1995) associated this 
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crossover with the transformation from the z — 1 quantum-critical to the z — 1 quantum- 
disordered regime. We have associated the same crossover with the development of the 
precursors of the spin- density- wave state in the electronic structure. We now show that the 
two identifications of the crossover at T* are in fact complementary. We note first that in 
order to obtain predominantly z — 1 behavior for the observables in the quantum-disordered 
regime, one must get rid of the Landau damping term in the dynamical staggered suscep- 
tibility. Otherwise, the system at the lowest T will necessarily crossover to the quantum- 
disordered z = 2 regime, where, e.g., 63 TiT oc uo s f = const, which does not agree with the 
data below T* (we recall that the measured m T\T increases as the temperature is lowered 
below T*). What happens when the quasiparticle spectrum develops a pseudogap near the 
corners of the Brillouin zone? At first glance, one might argue that the spin damping must 
become exponentially small in T at low temperatures because the former hot spots which 
chiefly contributed to the damping in the overdoped regime disappear. However, this argu- 
ment is incorrect because in a process of the Fermi surface evolution towards small pockets, 
there appear new hot spots which survive even in a situation in which the spin-density-wave 
structure of the electronic states is already developed. 

The actual reason the damping term goes down sharply when pockets are formed relates 
to the form of the fully renormalized spin-fermion vertex: in the spin- density- wave state with 
long-range order this vertex vanishes identically for the bosonic momentum Q because of 
the Ward identity (Adler principle); thus, the damping term in the spin susceptibility in fact 
scales as Jmx -1 (g, u ) °^ iui(q — Q) 2 (Sachdev 1994, Schrieffer 1995, Sachdev, Chubukov and 
Sokol 1995). In the precursor to the spin- density- wave state, there is no precise requirement 
that the full vertex at Q should vanish; however explicit calculations show that vertex 
corrections (which in the g ^> C limit have a negative sign) almost completely cancel the 
bare interaction g such that the fully renormalized vertex turns out to be small compared 
to g by a factor (g£)/g) 2 (Chubukov, Morr and Shakhnovich 1996). Very similar arguments 
have been previously provided by Schrieffer (Schrieffer 1995). 

A nice feature of the large g limit in the spin-fermion model is that the smallness of g^) / g 
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not only yields a small vertex but also allows one to compute the full vertex and full fermionic 
Green's function by expanding in powers of g'f) jg. This in turn allows one to compute 
fermionic damping: one just substitutes the full Green's functions and full vertex into the 
particle-hole bubble. We have performed this calculation and found that the damping term 
contains a small factor (g@) /g) A . Since the ratio gjgy decreases with T, the damping term 
rapidly decreases with decreasing temperature, and at low temperatures becomes smaller 
than CswC," 1 which is the energy scale of a gap in the spin susceptibility. In this situation, 
the system should display predominantly z = 1 quantum-disordered behavior, in agreement 
with the data below T*. In particular, u a f oc I/7 should increase as T decreases which in 
turn leads to the increase of e3 TiT with decreasing temperature. In addition, below T*, the 
resistivity scales as p oc T 2 /u s f, i.e., it rapidly (faster than T 2 ) decreases with decreasing T. 
This rapid decrease of p is also consistent with the experimentally observed behavior of the 
resistivity (see Fig|6])). Finally, at even lower temperatures, the system should, in principle, 
undergo a crossover to z = 2 quantum-disordered regime, but the crossover temperature is 
likely lower than T c which implies that this crossover cannot be observed experimentally. 

The reduction of the spin-fermion interaction vertex also affects the low T behavior of the 
uniform susceptibility. In a fully developed spin-density-wave state with long-range order 
only the longitudinal susceptibility has a (doping dependent) Pauli contribution, associated 
with a finite density of holes; the transverse susceptibility x± & t finite doping remains 
virtually the same as at half-filling because the Pauli contribution to x± is reduced by vertex 
renormalization (Chubukov and Frenkel 1992, Chubukov and Musaelian 1995). Specifically, 
we have 

X { I=° ] = ^ + 0(x); X? z = 0) = \x P (17) 

where x is the doping concentration, Z x is the quantum renormalization factor, and for a 
Fermi surface consisting of small pockets near (tt/2, vt/2), one finds Xp — y/mim2/2iT, where 
mi and irt2 are the two effective masses for fermionic dispersion. In order to estimate Xp-, 
we use the photoemission data for the oxychloride Sr 2 Cu0 2 Cl2 (Wells et al 1995, LaRosa 
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et al 1996). These data yields mi ss m 2 ~ 1/2 J, or %p ~ l/47rJ. In the preformed spin- 
density-wave state the transverse and longitudinal susceptibilities are indistinguishable so 
that 

11 Z 1 

X„(T -> 0) = -xl + -X^ = j^j + gXp- (18) 

Consider further the underdoped cuprates with a magnetically disordered ground state 
and a small Fermi surface. As the correlation length saturates to a finite value as T — > 0, 
the magnetic part of the susceptibility vanishes (Z x — > 0). We are thus left with only the 
Pauli contribution, i.e., Xu = (l/3)(Xp/2) ~ l/247rJ. We emphasize that the factor 1/3 in 
Xu is due to the spin-fermion vertex reduction, i.e., the same effect which yields a reduction 
in the spin damping. Reinserting the (g^B) 2 factor in Xu we obtain Xu/a*b ~ 0.4states/eV 
per Cu atom which is some six times smaller than Xul H 2 b — 2.6states/eV, the result for an 
optimally doped sample such as YBa 2 Cu 3 07, or for overdoped samples. We see therefore 
that the uniform susceptibility in the preformed spin-density-wave state with a small Fermi 
surface is much smaller than that found at larger doping levels with large Fermi surfaces, 
despite the fact that the Pauli susceptibility in 2D does not depend on p F and is the same 
for large and small elliptical Fermi surfaces. 



V. OPTIMAL DOPING 

Optimally doped cuprates are frequently defined as those which, within a given family, 
exhibit the highest T c . Examples of optimally doped materials include, e.g., YBa 2 Cu 3 06.93 
and Lai.85Sr 15 Cu0 4 . Photoemission experiments on YBCO systems and the bismuthates 
have shown that these materials have a large Fermi surface in close analogy with overdoped 
cuprates. At the same time, from a magnetic perspective, optimally doped materials are 
members of the underdoped family, in that their uniform susceptiblity exhibits the same 
crossovers at T cr and T* as are seen in other underdoped materials. The analogy between 
underdoped and optimally doped materials has also been observed in photoemission exper- 
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iments on Bi 2212 (Ding et al 1996): in both types of materials, there exists a pseudogap in 
the normal state which disappears at T ~ T cr which is substantially larger than T c . 

Another feature of optimally doped materials is that their resistivity continues to be 
linear in T down to T ~ T c . According to the arguments we have given above, such linearity 
above T reflects both the closeness of the Fermi surface to the magnetic Brillouin zone 
boundary and the particular spin-fluctuation spectrum. It is tempting to conjecture that in 
the optimally doped cuprates the Fermi surface is rather flat in such a way that (Ak max ) 2 
is minimized, which in turn yields a minimum in To as a function of doping. This flattening 
of the Fermi surface would also give rise to stronger vertex corrections which, as we recall, 
act to increase the spin-fermion interaction vertex and hence T c , before any changes in the 
Fermi surface topology take place. 

VI. SUMMARY 

We have proposed a specific scenario for the temperature crossovers in the overdoped 
and underdoped cuprates. We considered a nearly antiferromagnetic Fermi-liquid model 
and argued that in the overdoped regime the spin damping, 7 oc (cj s j£ 2 ) -1 is independent of 
temperature over the entire experimentally probed T-range. In this situation the magnetic 
behavior is described by a mean-field z = 2 dynamical susceptibility. We argued that the 
resistivity has a crossover from a Fermi-liquid like T 2 behavior at T < T to a linear in 
T behavior for T > T ; the crossover temperature T is low due to the proximity of the 
Fermi surface to the magnetic Brillouin zone boundary. We further argued that in order 
to account for the experimentally measured sequence of crossovers in underdoped cuprates 
from the Fermi-liquid, z = 2 behavior at high temperatures to the z = 1 scaling behavior 
at intermediate T, and to the pseudogap behavior at even lower T, one has to take into 
consideration the thermal variation of the damping rate of spin excitations, 7. We argued 
that the primary source of the variation of 7 with T in the underdoped cuprates is the thermal 
evolution of the quasiparticle dispersion near the Fermi surface produced by the precursor of 
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a spin-density-wave state. Within our approach, we found consistency with the experimental 
data on resistivity, NMR relaxation rates and uniform susceptibility in underdoped cuprates 
(for a discussion of the application of NAFLM to the Hall effect measurements see Stojkovic 
1996). We emphasize, however, that here we only discuss a possible scenario of how the 
damping evolves with T. At present, we can compute this damping either for very small 
or very large values of the coupling constant. Direct calculations of j(T) within the spin- 
fermion model at intermediate couplings are clearly called for. 

A final note. Our scenario is in contradiction with the proposals (Emery and Kivelson 
1995, Randeria et al 1994, Ding et al 1996) that the physical origin of the pseudogap behav- 
ior is the precursor to the d— wave pairing state. Both this and our scenario imply that the 
quasiparticle gap near (0, n) observed in photoemission measurements in the normal state of 
underdoped cuprates should not change as the system becomes superconducting, in agree- 
ment with the data. In the precursor to the d— wave pairing scenario, the superconducting 
gap is already preformed by superconducting fluctuations, while in our scenario there exists 
a preformed spin-density- wave gap near (0,7r). 

The photoemission data of Ding et al (1996) have been interpreted as an evidence in 
favor of preformed d-wave gap. They measured a gap as a shift of the midpoint of the 
leading edge of the photoemission spectrum, and found that the shift, which is of the order 
of 20meV, has the same /c-dependence as the d-wave gap in the superconducting state. On 
the contrary, Marshall et al (1996) estimated the gap from the position of the maximum 
in the spectral function. They found a much larger gap of about 0.2eV for 60K super- 
conductor, which is fully consistent with spin-density-wave scenario (in a fully developed 
spin-density-wave state, the gap near (0, tt) is about 2 J ~ 0.25eV). In the absence of a the- 
ory for the lineshape, it is difficult to say with confidence which interpretation of the data is 
correct. We note, however, that no preformed gap has been observed in overdoped cuprates. 
While the magnetic scenario leads naturally to the prediction that there should be no sign 
of a pseudogap in overdoped cuprates (superconductivity occurs before spin-density-wave 
precursor can develop), it is not clear why precursors to the d— wave pairing state should 
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disappear in overdoped materials. Finally, we note that the magnetic scenario correctly 
describes the whole sequence of crossovers in the normal state including the crossover at T cr 
which can easily be much larger than T c (e.g., T cr m 6T C in YE^CuiOs). Moreover, NMR 
experiments on YBa2Cus07, YE^CuztOg and Lai.ssSro.isCuO^ suggest that the magnetic 
correlation length is appoximately the same at T cr , £(T cr ) ~ 2a, in all of these compounds, 
implying yet another connection between the magnetic and the pseudogap behavior. It also 
seems unlikely that one could explain the crossover to z — 1 scaling as due to a precursor 
of the d-wave pairing. We therefore believe that the pseudogap has a magnetic rather than 
superconducting origin. 

Several researchers (see, e.g., Millis and Monien 1994) have argued that the pseudogap 
behavior in the underdoped cuprates is due to an exchange coupling, J±, between the bilayers 
found in a unit cell. Indeed, strong bilayer coupling leads to singlet configurations of adjacent 
spins which gives rise to a gap in the spin excitation spectrum. In an insulator, however, 
one needs a rather large Jj_ (~ 2. J) to produce the single configuration between adjacent 
spins in the bilayers, while recent experiments find Jj_ ~ 0.1J (Keimer et al 1996). For 
doped materials, Ioffe et al (1994) argued that the effective bilayer coupling scales as J e /* oc 
J±X 2 (Q), i-e., it is enhanced if susceptibility is strongly peaked at Q. This effect is certainly 
present at optimal doping, however we argued above that in the underdoped cuprates (where 
pseudogap has been observed), the enhancement of susceptibility is compensated by vertex 
corrections. In this sutuation, J e ^ should be of the same order as the bare coupling and is 
unlikely to give rise to a spin gap unless one assumes that there is a spin-charge separation 
(Millis, Altshuler and Ioffe 1996, Ubbens and Lee 1994). Our point of view is that the 
pseudogap behavior is an intrinsic property of a single Cu02 layer as the data on the spin- 
lattice relaxation rate, the uniform susceptibility, and £(T), in the single layer La 2 - 2: Sr x Cu04 
materials display the same sequence of crossovers seen in the bilayer YBa 2 Cu30 6+x materials. 
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FIGURES 

FIG. 1. The schematic behavior of various observables in underdoped cuprates (after Barzykin 
and Pines 1995). The temperature scales T cr and T* indicate crossovers which are defined and 
discussed in the text. 

FIG. 2. (a) Copper spin-lattice relaxation rate 63 T\T as a function of temperature in slightly 
overdoped Tl2Ba2Cu06+a; (Itoh et al 1994). Above 150K 63 TiT and, hence, uj s f is proportional to 
T. At lower temperatures, 63 TiT becomes T independent as in a conventional Fermi-liquid; (b) 
the ratio 63 TiT/T 2 2 G , where T^g is the spin-echo decay time. This ratio is practically independent 
of temperature, which is consistent with z = 2 scaling. 

FIG. 3. A sequence of crossovers in overdoped (a) and underdoped (b) cuprates. 

FIG. 4. The resistivity as a function of temperature for two overdoped Ti2Ba2Cu06+<5 samples 
(Kitaoka et al 1991) with T c = 40K (upper curve) and T c = OK (lower curve). The second 
sample is more heavily overdoped. Notice a pronounced T 2 behavior in both samples at very low 
temperatures. The top curve clearly exhibits a crossover to p oc T somewhere between 100 and 
200K. The crossover temperature for the bottom curve is presumably larger than 300i"C. 

FIG. 5. (a) The ratio 63 TiT/r| G as a function of temperature in underdoped YBa2Cu40s. As 
in overdoped samples, this quantity becomes independent of temperature above T cr ~ 450 — 500K; 
(b) 63 TiT/T2G as a function of temperature for the same compound. This ratio is almost a constant 
in the temperature interval T* < T < T cr , which is a signature of the z = 1 scaling regime. In both 
graphs, the dashed line is a guide to the eye. 

FIG. 6. The resistivity along the a— axis for underdoped YBa2Cu40s- Notice that the linear 
behavior starts at around T*, and the minor change of slope at T = T cr . Below T*, the resistivity 
rapidly falls off with decreasing T. 
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